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^ ; I. INTRODUCTION 

O; 

The problem of strongly correlated electrons in one dimension attracted the attention of condensed matter physicists 
(N . since early 50's. After the seminal paper of Tomonaga Jl|] who suggested the first soluble model in 1950, there were 
essential achievements reflected in papers by Luttinger 0] (1963) and Mattis and Lieb || (1963). Considerable 
contribution to understanding of the generic properties of the ID electron liquid has been done by Dzyaloshinskii and 
Larkin [Q (1973) and Efetov and Larkin H (1976). Despite the fact that the main features of the ID electron liquid 
. were captured and described by the end of 70's, the investigators felt dissatisfied with the rigour of the theoretical 
description. The most famous example is the paper by Haldane || (1981) where the author developed the fundamentals 
, of the modern bosonisation technique, known as the operator approach. This paper became famous because the author 
D ■ has rigourously shown how to construct the Fermi creation/anihilation operators out of Bose ones. The most recent 
^ \ example of dissatisfaction is the review due to von Delft and Schoeller (1998) who revised the approach to the 
■ bosonisation and came up with what they called constructive bosonisation. 

There exists another way to bosonise the problem. It was suggested 25 years ago by || and found its further 
elaboration in the paper by Lee and Chen || in 1988. This approach (sometimes called "functional bosonisation") 
can be found in variety of papers (see, for example, |l0| ] and JTT|). The main objective of this lecture is to assemble 
these pieces in a clear picture. The aim is purely methodological which is the reason why the number of complications 
""O ■ (like backscattering, presence of lattice and spin degrees of freedom) are kept at the lowest possible level. Any 
generalisations can be easily incorporated into the scheme. The final formulation of the bosonised problem will look 
distinct from what can be found elsewhere. We believe that the approach present in this lecture is the "shortest 
shortcut" from the original fermionic problem of interacting electrons to the problem of non- interacting bosons. 
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II. FUNCTIONAL INTEGRAL FORMULATION 
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We start from a general Hamiltonian describing interacting electrons in one dimension: 



ii = J— [ dxd x &(x)d x $(x) + I [ dxdx'&(x)&(x') V (x - x') $(x')$(x). (1) 
2m J 2 J 

Here Vo(x — x') is the bare inter-electron interaction which is usually taken to be the Coulomb interaction screened 
by the carriers of the media surrounding the one-dimensional system. 

The thermal Green function is the time-ordered product of the Fermion anihilation and creation operators 

a! g(x,T;x',T') = {T T *{x,T)*(x',T')) (2) 

o 
o 



defined in the imaginary time Heisenberg representation 

$(x, t) = e^' r *(x)e^' r , V(x, r) = e^' T ¥ {x)e^' T , (3) 

where TL' = TL — fiN. Here /i is the chemical potential and N is the particle number operator. This Green function can 
be written in the form of the coherent state functional integral over the Grassmann fields that reflects the fermionic 
nature of electrons: 

y V> T > x > T ) /XtyX>exp{-ShM*]} ' {) 

where the action S[ip, ip*] is defined by 

5[V>,V>*]= / dxdT^*(x,r) d T +i ip{x,T) + ]- I dxdx'dTi/;*(x,T)Tp*(x',T)Vo(x-x')Tp(x',T)ip(x,T). (5) 
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Here £ — — (2m)~ 1 d^, — fi is the kinetic energy operator counted from the chemical potential /i, and ip(x : r) is the 
Grassmann function (or, equivalently, fermionic field) which is an anti-periodic function of the imaginary time r, 
i/j(t + 0) = —i/j(t), with the inverse temperature (3 = T _1 being the period. 

Being aware that individual electrons are not good as a starting point for description of excitations in one dimension 
we choose to introduce an auxiliary field conjugated to the electron density. This can be achieved via the Hubbard- 
Stratonovich transformation which brings the fermionic problem to a Gaussian form at the price of complication 
caused by appearance of some new field and related to it extra integration. We start from the following representation 
of the exponential containing the interaction 



exp [Sint] 
where we used symbolic expressions 



exp [-i^fV] 



(6) 



0VV= / dadx'dr^(a:,T)V ~ 1 (a!-a/)^(a/,T), 



(jnp*ip = J dxdr(j}{x, r)t/j* (x, r)ip(x, t) 

The notation V^~ 1 (x — x') stands for the inverse function understood as an integral kernel inverse to that corresponding 
to the bare interaction: 



dxiVjj (x — xi)Vq(xi — x') — S(x — x'). 



(7) 



One can see that the field (f>(x, r) should posess the same properties as a product of two Grassmann fields 
ip*(x,T)ip[x,T), i.e. it is a commuting field and periodic under the shift r — > r + f3. Therefore, this field is of 
bosonic nature. Now we substitute the Hubbard- Stratonovich transform, Eq.(^J), into the original representation 
Eq.([|) of the thermal fermionic Green function: 



G(x,t;x',t') 



J V^Vij)*V(j)^(x, r)ip*(x', t') exp {-S[ip, 0]} 



fVtl)T>i/}*T>^exp{-S[ilj, ip*;c/>\} 
Now the action is a functional of not only the fermionic field but also of the auxiliary field 0: 

S[ip , ip* ; (/>] = — J dxdx'dT<j)(x,T)V ( f 1 (x — x / )(f>(x / ,T)+ J dxdr?p*(x, r) d T + £ — i(f>(x, r) tp(x,r) 

Now we divide and multiply both numerator and denominator by 



DipD'ip* exp 



dxdr-0* 



d T + £ - 1$ 



V> } = Det 



d T + £ - ic. 



— exp VTr In 



d T + £ - iq 



(8) 
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Then one can see that the thermal Green function can be written as some auxiliary Green function G[4>] 
g(x,r;x',r';[ct>]) 



J VipVip*^{x,T)ip* (x' ,t') exp ^- J dxdrip* 




\-i-i<t> 




J VipV-ip* exp j — J dxdrip* 


dr+i- 


- icf> 





averaged over the field <p with the new bosonic action S[4>]: 

g{x,T-,x' ,t') = z~ l j v^g{x, 



t; x',t'; 



where 



Z = / V6e 



S[4>] = ~ Tr ln °r + i - i<t 



(11) 



(12) 



(13) 



The auxiliary Green function is the one of the "free" electrons in the spatial and temporal dependent field <p(x,r) 
and it must be found from the equation 
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(d T +i-i4>(x,T^j g(x,T;x',T';[<p])=5(x-x')5(T-T'). 



(14) 



If we were able to solve this equation for arbitrary field 4>(x, r) we would then have some expression for Q which 
depends only on the bosonic field </>(x, r). Substituting this solution of Eq.(|l4|) into Eq.(p"^) we would finally reduce 
our problem to another one where we must do some integral over the bosonic field governed by the action Eq. (fl3|) . 
And it should be pointed out that this procedure can be implemented for any observables constructed from any 
number of Fermi operators, not only for a bilinear combination defining the one-particle Green function Eq.(^). The 
only extra effort we would need is to use the Wick theorem in Eq.(|li"|) which is valid because the action is quadratic 
in V'-ficlds after the Hubbard-Stratonovich transformation is implemented. 

It is evident enough, that this scheme is absolutely general and therefore absolutely useless. We should refer to 
some specifics of one-dimensional problems to develop it to the stage where it becomes a working tool. And that is 
what we will do next. 

III. ONE DIMENSION 

The main feature of one dimensional problems with Fermions is that their Fermi surface consists of just two points 
— pf and +pf- Accordingly, electrons close to this Fermi surface can be divided into two groups, left- (with the 
momenta in the vicinity of —pf) and right-movers (with the momenta in the vicinity of —p F ). That means that the 
Fermi fields over which we integrate in Eq.(|]) might be represented in the form 

4>(x,t) =^+(x,T)e tpFX + ^-(x,T)e- tpFX , (15) 

We are interested in large scale behaviour of our system which makes this distinction well defined because only small 
vicinity of the Fermi points should be taken into account. In other words, the Fermi fields ip±(x, r) may be considered 
as smooth functions on the scale Pp 1 ■ It allows as when writing action in Eq.(Q) to act by the kinetic energy operator 
£ on the fast oscillating exponential mainly, neglecting second order derivatives of i/>±{x,t). Also, the integral over 
the whole space suppresses contribution from the "uncompensated" oscillatory functions proportional to e ±l2pii,a: . 
Therefore, we may introduce two-vector 

V> t {x,t) = (i/> + (x,t) 1 i/>-(x,t)) (16) 
and in these notations the action Eq.(|l^) becomes 

i Tviv, ( 9 T - iv F d x ~ i(j> 



Writing down this action we neglected the backscattering of electrons due to Coulomb interaction. Generally speaking 
in the Coulomb field <f)(x, r) there are two essential components related to two different scattering events. One of 
them contains the Fourier harmonics with \q\ <C Pf and, therefore, is a smooth function of x, another one consists of 
the Fourier harmonics with |q| « 2pp- The former is responsible for the forward scattering, i.e. when two colliding 
electrons keep their direction of motion which they had before collision. The latter causes backscattering when two 
colliding electrons make the U-turn after the collision. If we assume that for the Fourier transforms of the bare 
interaction we have the inequality Vo(q <C %Pf) *C Vo(q ~ 2pp ), we may neglect the backscattering which means that 
the Coulomb field <j>{x, t) in the action Eq.(|l7|) is a smooth function on the scale of the electron wavelength pp 1 . 

The Green function now becomes a 2 x 2 matrix and neglecting backscattering we find that it is a diagonal matrix 
with the diagonal elements Q+ and (?_ defined by the equations: 

{d T T iv F d x - i<f>(x, t)) G ± (x, t; x', t'; [<£]) = S(x - x 1 ) S(t - t'). (18) 

The formal solution of these equations can be easily written 

g + (x, r; x',t'; [0]) = g + (x -x',r- r') e ^(^M^V). (19) 
r; x',t'; [</>]) = 9 -(x - x', r - r') e *>*(*>T)-*> W), 

where 9(x, r) should be determined from another equation 

(d T -iv F d x )6{x,T)=cl>(x,T) (20) 
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and g± = (d T ± iv F d x ) 1 are the Green functions of the chiral particles: 

9±{x,t) = JL S m- l KT[TTi— )■ (21) 
2v F \ v F J 

Now, when we have conveniently expressed the auxiliary Green function as an exponential of the linear functional 
of the 0- field, one should take a closer look at the action itself. The term written as Trln in Eq.(^) is a sum of two 
terms which are the contributions of the right- and left-movers. Let us analyze the right-movers contribution first. 
We may expand Tr In in powers of <fi: 

\j * 

Trln (d T - iv F d x - Mx, r)) = Trln (d T - iv F d x ) - V -Tr (g+ icj)) n . (22) 

n—l 

The term containing the n-th power of <j> is the vertex containing the loop of n-th order F+ and n external lines 
corresponding to </>'s: 

r. n n 

Tl '(5+0)"= / JJ dx k dr k r+(xi,n; ...;x n ,T n )H<Kx i ,T i ). (23) 

fc=l i=l 

The loop is built from chiral Green functions g + in the following way: 

n 

r^(xi,ri;...;a: n ,T n ) = - x i+ i,Ti - t 1+1 ) (24) 

i=l 

and (x„+i, t„+i) = (xi,Ti). It is convenient to calculate this loop in the frequency/momentum representation. We 
use notations e and fl for the fermionic and bosonic Matsubara frequencies correspondingly. Due to translational 
invariance of the system the loop T„ will depend on n — 1 momenta and frequencies 

f dp 

Using explicit expression of the Green function of right-moving free electrons 

9+(p, e) = + v F q)~ (26) 

and introducing temporarily the notations 

z = — ie + v F q, z% — — ifij + v F qt (27) 
we can write the expression for the n-th loop in the following form 

r£ = / dz — r, (28) 

7 Z (z + 21 ) (z + Zl + Z2) — {Z + Zl+ Z 2 + ... + Z n -l) 

where we have used a symbolic notation for J dz — T^] £ J There are very widely used in QED the Feynman 
identities which help calculating higher order diagrams. One of them, 

I _ f (n- l)!dxi...dx n _i 



ai...a 2 J [aix n ~i + a 2 (x n ^ 2 ~ x n -i) + ... + a n (l - xi)] n 

0<a„_i<...<xi<l 

after some modification, may be brought to the form which is of help in our case, namely 

n-1 



r+(zi,...,z„_i) = (n-l)! / dz / d n - L y 



z + ^ViZi 



i=l 



(30) 



where the inner integral runs over the area Y n -\ = (0 < y n -i < ••■ < j/i < 1). All loops r+ appear only in integrals 
with the products of </>'s (see Eq.(p3[)) which are explicitely symmetric with respect to the permutation of all argu- 
ments (xi,Ti) (or, equivalently, with respect to the permutations of in frequency /momentum representation). 
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Therefore, only symmetrized form Tn ^ sym ^ contributes to the integral in Eq.d23l). Symmetrization of of the expression 



( |30| ) leads to discarding the factorial and relaxing the ordering of j/'s: 



r+ (sym) (*i,-,^-i) = Jdz J d* 



n-l 



1=1 



(31) 



For n > 2 these integrals vanish. To see that we may take one integral over some j/i, say y n -V 



Y+ (syn 



dz I d n - 2 y 



Z + Zr, 



n-2 



-n+1 



n-2 

E 



(32) 



For n > 2 the integral of each term in the square brackets is finite. It means that we can make a shift of the integration 
variables in each term independently. Making the shift z + z„_i — > z (meaning p + q n -i — > p, e + fi„_i — ► e) in the first 
term (integral) we find that it becomes identical to the second one and, thus, cancel it. Therefore, all loops containing 
more than two external lines are zero. That was the observation first made by Dzyaloshinskii and Larkin [W. Above 
I gave the simplest proof to that fact. Nowdays this cancellation when it happens in many different one-dimensional 
problems is called the Loop Cancellation theorem. 

Therefore, we are left with only two contributions from the loops with one and two legs. Calculating them one 
must excercise precaution because the corresponding integrals are ill-defined. These troubles have their roots in the 
linearization of the spectrum. However, we may predict the correct result even without healing these problems. It 
is clear that the loop with one leg is proportional to the zero mode of the Coulomb interaction and unperturbed 
particle density. This term is always cancelled against positive background. The loop with two legs is the standard 
polarization operator which is know to be (for the right-movers) in (g, ^-representation 



1 



2tt —id + vfQ 



(33) 



Finally, back to the action Eq.(|17|), we collect two contributions coming from both chiral species and notice that the 
part independent of <f> drops out of the expression Eq. (fl4|) . The action then becomes 



(34) 



where we introduced the overall polarization operator which is the sum of two chiral polarization operators 7r_|_ and 
7r_ and has the following familiar form in (g, ^-representation 



n(q, ft) = TT+(q, fi) + ?r_(g, fi) = v 



2 2 

v F q 



Q 2 + v 2 F q 2 ' 



(35) 



with v = (ttvf) -1 being the one-dimensional density of states. 

The Gaussian action describing Coulomb bosons <f> can be cast in a familiar form if we introduce effective 
time / frequency-dependent interaction 



s[<t>] ^^v- 1 ^ 



(36) 



where 



V~ 



Vn- 1 + n, 



V- 1 (q,n)=V - 1 (q)+K(q,n) 



(37) 



One can easily recognize the well-known result for the effective Coulomb potential with the screening calculated in 
the random phase approximation. It means that the random phase approximation becomes exact as we linearize the 
spectrum. 



IV. CALCULATION OF OBSERVABLES 

After the Hubbard-Stratonovich transformation is done the action Eq.(^J) becomes quadratic in the fermionic fields 
and Wick's theorem becomes applicable for any product of the ^-fields. Consequently, any observable, originally a 
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product of ^-operators, decomposes into the product of the auxiliary Green functions Q± [(/)] which should be averaged 
over the auxiliary field <f>. This field describes non-interacting bosons whose propagator is the screened Coulomb 
interaction Eq. ( pq ) . The auxilary Green functions Q± [</>] are products of the free electron Green function g± and the 
exponentials with the exponent being a linear functional of the Coulomb bosons (j> (see Eqs.(|l9|)and (po|)). Therefore, 
the problem of interacting electrons is reduced to the calculation of the simple Gaussian integrals. 
In this way we, as an example, may calculate the Green function of the right movers 



where the exponent is given by 



G+(x,t; [</>])) = .g+(a;,T),exp [-K(x,t)\ 



dq V(q,Si) 
2tt (-ifi + v F q) 2 



1C(x,t) = 1{{9(x,t)-6(0,0)} 2 )=tJ2 I 



iqx—iQr 



The propagator of the "screened Coulomb field" (0-field) is 

n 2 + v 2 F q 2 



V(q,Q)=V (q) 



n 2 + v 2 q 2 ' 



where the renormalized velocity v is defined as 



v(q) = v F \ 1 + 



Vo(q) 
ttvf 



The simple arithmetics gives us the result which can be found in any review on one-dimensional bosonisation: 



K(x, t) 



dq 
<1 



e —v F qr 



-vqr 



e -0v F q _|_ I 



-fivq _|_ I 



iqx 



(v — vp) 2 1 — cos qxe vqT 



2vvf 



-fJvq _|_ I 



(38) 



(39) 



(40) 



(41) 



(42) 



It is not our goal to describe the behaviour of different response functions. The results obtained with the use of the 
formalism present in this lecture are exactly the same as those obtained within operator formalism due to Haldane. 
The main distinction of the approach given in this lecture is its simplicity. 



V. ELECTRONS WITH SPIN 



We accept "g-ology" as our starting point, neglecting Umklapp and backscattering: 

S[tp] = / dxdTip* a (x, t) [d T - ir a v F d x ] VVa^, r) 



+- I dxdx'drg^, {x - x')n aa {x,T)n a i a >{x' ,t), 



(43) 
(44) 



where a —\ , J, describes spin dependence and a = ± refers to right- (+) or left-moving (— ) particles (r± = ±1). The 
coupling matrix has the structure 



Sera 1 — 3a,a' $<T,<j'g\ + 5a,-a'9i 

After Hubbard-Stratonovich transformation we arrive at 

s i4>] = \<t>s9^(t>s' - y^Trln [d T - ir a v F d x - i(j> s ] 



(45) 



(46) 



Please note that we combine spin and chirality in one composite index s = (c, a), where possible. Then expanding 
'Trln' as in the previous section we are again left with the RPA result: 



S[4>] = 5rpa[</>] = 7^<t>sV ss }<p s >, V ss } = g ss ) + S s>s ,TT a 



(47) 



G 



The matrix Green function is diagonal in the spin indices and is the sum of two chiral functions: 

Sao'Cz-a/.T-T'J^^^afoT)^ (48) 
a, a' a 

where 

Q aa {x -x',t- t') = g a {x - x',r- t') ( e iV>,(*,T)-e s (x>,T>)]^ (4Q) 

The brackets stand for the averaging over 9 S with the action Eq.(f47|). The field 9 S itself is a linear functional of the 
"Coulomb field" <j) s , and must be found as the solution to the following equation 

(d T - ir a v F d x ) S ( X ) = (p s (x). (50) 

Taking the Gaussian integral over <j> we obtain 

' e i[0 3 {x,T)-6 3 {x' ,t')]\ _ e -X;,(x-x',r-T') j ( 51 ) 

with the exponent related to the diagonal matrix element of the screened 'Coulomb interaction' V BS (q, fi) (see Eq. ([fTf) ) : 

Due to the fact that the free Hamiltonian (and thus, density-density correlators 7r±(q)) do not depend on spin, we 
can make a rotation in the spin sub-space only to diagonalize the coupling matrix g ss i . The rotation can be performed 
by the matrix which has the form 

L = ^(!-i) (53) 

in the representation where the four-vector composed of components <f> ua is taken to be (j) T — (<p^+, <fi-\-, <Pi+i 4*1— )• 
This rotation, in fact, is equivalent to saying that we are writing the quadratic form in the new fields <p pa and <p va 
which arc the linear combinations of the old ones: 



-^=(0 Tq + iQ ), <p aa = 



' J P a - ~7^(^ a + 4 la), ^aa = "a^T" _ 'Plot), (54) 



corresponding to charge- (index p) and spin- (index a) fields. After such a rotation the action becomes 

S[cp] = \ Y, 4>^V-^ a >, (55) 

v, a, a' 

where v = p,a. The new interaction matrix is block-diagonal in its spin/charge sectors and each block can be 
represented as a 2-by-2 matrix in the chiral indices: 

% = - ( , fa 4u + ^ a 2 \ 

{a iv + K + ){ai v + 7T_) - a|„ \ a 2v a 4v + n+ J 

where we introduced spin and charge coupling constants 



9i P = \{g\ + 9i~)i 9i* = lid! ~ 9^), a tu = of 2 9w _ 2 y (57) 

z z z \9iu 9lv) 



Substituting explicit expressions for the chiral polarization operators ir±(q, Q) one can write it in the following form: 

MqM) = 2n ^ 2 + V l q2 2 A(q,n) (58) 
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where A(q, Q) is 2-by-2 matrix 

A(q,n) = 



( -iQ + v F q 
1 (y - v„) - 2~f„v v 



V 



ifl + vpq 

-%{K V - K^) (v v v ) 2 lv v v — f + r '- q 



—iQ + vpq 



written in the notations standard for the literature on the Luttinger Liquid: 



vt = Vf + ~{gi v ± 92u), K v = 



(59) 



v v = V vtvu , 7„ = - (K v + K v 1 - 2) 



The diagonal part of the interaction matrix V ss needed for the average Green function does not actually depend on 
spin index a: 



Therefore, for the average Green function we have the representation: 

Q{x,t) =Y,9 a {x,T)e- K ^e i ^ x , K+{x,t) = K-(-x, -r) 



(60) 



(61) 



where 



\ _ e -v„qr+iqx ^ 



-Vpqr+iqx 



e -f3v F q + 1 



e -vuqr ^| 

— 7i/ — 7a cos qx > 



(62) 



Once again, the expression above is well known and can be found in any review on the Luttinger Liquid. The task 
of this Lecture was to show that starting from the field-theoretical description of one-dimensional interacting electrons 
( by this I mean representation of the observables as functional integral over the Grassmann fields ) one can derive 
the free-bosons action in extremely simple, almost trivial, and elegant way. 
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